Cis-regulatory interfaces reveal the molecular mechanisms underlying the notochord gene regulatory network of Ciona.

Tissue-specific gene expression is fundamental in development and evolution, and is mediated by transcription factors and by the cis-regulatory regions (enhancers) that they control. Transcription factors and their respective tissue-specific enhancers are essential components of gene regulatory networks responsible for the development of tissues and organs. Although numerous transcription factors have been characterized from different organisms, the knowledge of the enhancers responsible for their tissue-specific expression remains fragmentary. Here we use Ciona to study the enhancers associated with ten transcription factors expressed in the notochord, an evolutionary hallmark of the chordate phylum. Our results illustrate how two evolutionarily conserved transcription factors, Brachyury and Foxa2, coordinate the deployment of other notochord transcription factors. The results of these detailed cis-regulatory analyses delineate a high-resolution view of the essential notochord gene regulatory network of Ciona, and provide a reference for studies of transcription factors, enhancers, and their roles in development, disease, and evolution.

Snail Transcriptionally Represses Brachyury to Promote the Mesenchymal-Epithelial Transition in Ascidian Notochord Cells.

Mesenchymal-epithelial transition (MET) is a widely spread and evolutionarily conserved process across species during development. In Ciona embryogenesis, the notochord cells undergo the transition from the non-polarized mesenchymal state into the polarized endothelial-like state to initiate the lumen formation between adjacent cells. Based on previously screened MET-related transcription factors by ATAC-seq and Smart-Seq of notochord cells, Ciona robusta Snail (Ci-Snail) was selected for its high-level expression during this period. Our current knockout results demonstrated that Ci-Snail was required for notochord cell MET. Importantly, overexpression of the transcription factor Brachyury in notochord cells resulted in a similar phenotype with failure of lumen formation and MET. More interestingly, expression of Ci-Snail in the notochord cells at the late tailbud stage could partially rescue the MET defect caused by Brachyury-overexpression. These results indicated an inverse relationship between Ci-Snail and Brachyury during notochord cell MET, which was verified by RT-qPCR analysis. Moreover, the overexpression of Ci-Snail could significantly inhibit the transcription of Brachyury, and the CUT&Tag-qPCR analysis demonstrated that Ci-Snail is directly bound to the upstream region of Brachyury. In summary, we revealed that Ci-Snail promoted the notochord cell MET and was essential for lumen formation via transcriptionally repressing Brachyury.

Conserved enhancers control notochord expression of vertebrate Brachyury.

The cell type-specific expression of key transcription factors is central to development and disease. Brachyury/T/TBXT is a major transcription factor for gastrulation, tailbud patterning, and notochord formation; however, how its expression is controlled in the mammalian notochord has remained elusive. Here, we identify the complement of notochord-specific enhancers in the mammalian Brachyury/T/TBXT gene. Using transgenic assays in zebrafish, axolotl, and mouse, we discover three conserved Brachyury-controlling notochord enhancers, T3, C, and I, in human, mouse, and marsupial genomes. Acting as Brachyury-responsive, auto-regulatory shadow enhancers, in cis deletion of all three enhancers in mouse abolishes Brachyury/T/Tbxt expression selectively in the notochord, causing specific trunk and neural tube defects without gastrulation or tailbud defects. The three Brachyury-driving notochord enhancers are conserved beyond mammals in the brachyury/tbxtb loci of fishes, dating their origin to the last common ancestor of jawed vertebrates. Our data define the vertebrate enhancers for Brachyury/T/TBXTB notochord expression through an auto-regulatory mechanism that conveys robustness and adaptability as ancient basis for axis development.

Genome-wide identification of notochord enhancers comprising the regulatory landscape of the brachyury locus in mouse.

The node and notochord are important signaling centers organizing the dorso-ventral patterning of cells arising from neuro-mesodermal progenitors forming the embryonic body anlage. Owing to the scarcity of notochord progenitors and notochord cells, a comprehensive identification of regulatory elements driving notochord-specific gene expression has been lacking. Here, we have used ATAC-seq analysis of FACS-purified notochord cells from Theiler stage 12-13 mouse embryos to identify 8921 putative notochord enhancers. In addition, we established a new model for generating notochord-like cells in culture, and found 3728 of these enhancers occupied by the essential notochord control factors brachyury (T) and/or Foxa2. We describe the regulatory landscape of the T locus, comprising ten putative enhancers occupied by these factors, and confirmed the regulatory activity of three of these elements. Moreover, we characterized seven new elements by knockout analysis in embryos and identified one new notochord enhancer, termed TNE2. TNE2 cooperates with TNE in the trunk notochord, and is essential for notochord differentiation in the tail. Our data reveal an essential role of Foxa2 in directing T-expressing cells towards the notochord lineage.

Eomes restricts Brachyury functions at the onset of mouse gastrulation.

Mammalian specification of mesoderm and definitive endoderm (DE) is instructed by the two related Tbx transcription factors (TFs) Eomesodermin (Eomes) and Brachyury sharing partially redundant functions. Gross differences in mutant embryonic phenotypes suggest specific functions of each TF. To date, the molecular details of separated lineage-specific gene regulation by Eomes and Brachyury remain poorly understood. Here, we combine mouse embryonic and stem-cell-based analyses to delineate the non-overlapping, lineage-specific transcriptional activities. On a genome-wide scale, binding of both TFs overlaps at promoters of target genes but shows specificity for distal enhancer regions that is conferred by differences in Tbx DNA-binding motifs. The unique binding to enhancer sites instructs the specification of anterior mesoderm (AM) and DE by Eomes and caudal mesoderm by Brachyury. Remarkably, EOMES antagonizes BRACHYURY gene regulatory functions in coexpressing cells during early gastrulation to ensure the proper sequence of early AM and DE lineage specification followed by posterior mesoderm derivatives.

Transcriptomic evidence for Brachyury expression in the caudal tip region of adult Ptychodera flava (Hemichordata).

Most metazoans have a single copy of the T-box transcription factor gene Brachyury. This gene is expressed in cells of the blastopore of late blastulae and the archenteron invagination region of gastrulae. It appears to be crucial for gastrulation and mesoderm differentiation of embryos. Although this expression pattern is shared by most deuterostomes, Brachyury expression has not been reported in adult stages. Here we show that Brachyury of an indirect developer, the hemichordate acorn worm Ptychodera flava, is expressed not only in embryonic cells, but also in cells of the caudal tip (anus) region of adults. This spatially restricted expression, shown by whole-mount in situ hybridization, was confirmed by Iso-Seq RNA sequencing and single-cell RNA-seq (scRNA-seq) analysis. Iso-Seq analysis showed that gene expression occurs only in the caudal region of adults, but not in anterior regions, including the stomochord. scRNA-seq analysis showed a cluster that contained Brachyury-expressing cells comprising epidermis- and mesoderm-related cells, but which is unlikely to be associated with the nervous system or muscle. Although further investigation is required to examine the roles of Brachyury in adults, this study provides important clues for extending studies on Brachyury expression involved in development of the most posterior region of deuterostomes.

The evolutionary history of Brachyury genes in Hydrozoa involves duplications, divergence, and neofunctionalization.

Brachyury, a member of T-box gene family, is widely known for its major role in mesoderm specification in bilaterians. It is also present in non-bilaterian metazoans, such as cnidarians, where it acts as a component of an axial patterning system. In this study, we present a phylogenetic analysis of Brachyury genes within phylum Cnidaria, investigate differential expression and address a functional framework of Brachyury paralogs in hydrozoan Dynamena pumila. Our analysis indicates two duplication events of Brachyury within the cnidarian lineage. The first duplication likely appeared in the medusozoan ancestor, resulting in two copies in medusozoans, while the second duplication arose in the hydrozoan ancestor, resulting in three copies in hydrozoans. Brachyury1 and 2 display a conservative expression pattern marking the oral pole of the body axis in D. pumila. On the contrary, Brachyury3 expression was detected in scattered presumably nerve cells of the D. pumila larva. Pharmacological modulations indicated that Brachyury3 is not under regulation of cWnt signaling in contrast to the other two Brachyury genes. Divergence in expression patterns and regulation suggest neofunctionalization of Brachyury3 in hydrozoans.

FHOD1 is upregulated in glioma cells and attenuates ferroptosis of glioma cells by targeting HSPB1 signaling.

BACKGROUND: As a new type of regulatory cell death, ferroptosis has been proven to be involved in cancer pathogenesis and therapeutic response. However, the detailed roles of ferroptosis or ferroptosis-associated genes in glioma remain to be clarified. METHODS: Here, we performed the TMT/iTRAQ-Based Quantitative Proteomic Approach to identify the differentially expressed proteins between glioma specimens and adjacent tissues. Kaplan-Meier survival was used to estimate the survival values. We also explored the regulatory roles of abnormally expressed formin homology 2 domain-containing protein 1 (FHOD1) in glioma ferroptosis sensitivity. RESULTS: In our study, FHOD1 was identified to be the most significantly upregulated protein in glioma tissues. Multiple glioma datasets revealed that the glioma patients with low FHOD1 expression displayed favorable survival time. Functional analysis proved that the knockdown of FHOD1 inhibited cell growth and improved the cellular sensitivity to ferroptosis in glioma cells T98G and U251. Mechanically, we found the up-regulation and hypomethylation of HSPB1, a negative regulator of ferroptosis, in glioma tissues. FHOD1 knockdown could enhance the ferroptosis sensitivity of glioma cells via up-regulating the methylated heat-shock protein B (HSPB1). Overexpression of HSPB1 significantly reversed FHOD1 knockdown-mediated ferroptosis. CONCLUSIONS: In summary, this study demonstrated that the FHOD1-HSPB1 axis exerts marked regulatory effects on ferroptosis, and might affect the prognosis and therapeutic response in glioma.

Uterine Sarcoma With FGFR1-TACC1 Gene Fusion: A Case Report and Review of the Literature.

With the growing availability of RNA sequencing technology in the pathology laboratory, new gene fusion-associated malignancies are increasingly being characterized. In this article, we describe the second ever reported case of a uterine sarcoma harboring a FGFR1-TACC1 gene fusion. The patient, a 53-yr-old perimenopausal woman, was found to have a 6 cm mass spanning the lower uterine segment and endocervix. Histologically, this was a spindle cell neoplasm with coagulative necrosis, moderate cytologic atypia, and increased mitotic activity. By immunohistochemistry, the neoplastic cells coexpressed CD34 and S100, and lacked smooth muscle marker expression. RNA sequencing revealed the presence of a FGFR1-TACC1 gene fusion. This report provides further evidence to suggest that FGFR1-TACC1 may be a recurrent fusion in a subset of uterine sarcomas. RNA sequencing using a panel that includes FGFR-TACC family fusions should be considered for uterine sarcomas that do not fit conventional diagnostic criteria, particularly as tumors with these fusions may be amenable to targeted therapy.

FGFR1/MAPK-directed brachyury activation drives PD-L1-mediated immune evasion to promote lung cancer progression.

Immune checkpoint inhibitors provide promising benefits for patients with cancer. However, efficacy has been encumbered by high resistance rates. It is critical to understand the basic mechanisms of tumor-mediated resistance to this treatment modality. Previous studies have found that the transcription factor brachyury is highly expressed in lung cancer. Here, we show that brachyury activation induces the upregulation of PD-L1 leading to inactivation of T cell proliferation in vitro and inhibited infiltration of CD8+ and CD3+ T cells into tumor in an immunocompetent mouse model. We further demonstrate that FGFR1/MAPK activation regulates brachyury and PD-L1 expressions and promotes immunosuppression. Blocking FGFR1/MAPK suppresses brachyury and PD-L1 expressions, revives immune activity, and reverses the resistance to anti-PD-1 treatment to produce a durable therapeutic response. We also find that lung cancer patients with high activation of the FGFR1-MAPK-brachyury-PD-L1 signature and low expression of CD8A, CD3D, or PDCD1 have worse survival outcomes. These findings elucidate a novel mechanism of immune escape from immune checkpoint therapy and provide an opportunity to enhance its therapeutic efficacy in the treatment of a subset of FGFR1/MAPK/brachyury/PD-L1-driven lung cancer.

Proteome profiling of enzalutamide-resistant cell lines and serum analysis identified ALCAM as marker of resistance in castration-resistant prostate cancer.

Enzalutamide (ENZA) is a frequently used therapy in metastatic castration-resistant prostate cancer (mCRPC). Baseline or acquired resistance to ENZA have been observed, but the molecular mechanisms of resistance are poorly understood. We aimed to identify proteins involved in ENZA resistance and to find therapy-predictive serum markers. We performed comparative proteome analyses on ENZA-sensitive parental (LAPC4, DuCaP) and -resistant prostate cancer cell lines (LAPC4-ENZA, DuCaP-ENZA) using liquid chromatography tandem mass spectrometry (LC-MS/MS). The top four most promising candidate markers were selected using bioinformatic approaches. Serum concentrations of selected markers (ALCAM, AGR2, NDRG1, IDH1) were measured in pretreatment samples of 72 ENZA-treated mCRPC patients using ELISA. In addition, ALCAM serum levels were measured in 101 Abiraterone (ABI) and 100 Docetaxel (DOC)-treated mCRPC patients' baseline samples. Results were correlated with clinical and follow-up data. The functional role of ALCAM in ENZA resistance was assessed in vitro using siRNA. Our proteome analyses revealed 731 significantly differentially abundant proteins between ENZA-sensitive and -resistant cells and our filtering methods identified four biomarker candidates. Serum analyses of these proteins revealed only ALCAM to be associated with poor patient survival. Furthermore, higher baseline ALCAM levels were associated with poor survival in ABI- but not in DOC-treated patients. In LAPC4-ENZA resistant cells, ALCAM silencing by siRNA knockdown resulted in significantly enhanced ENZA sensitivity. Our analyses revealed that ALCAM serum levels may help to identify ENZA- and ABI-resistant patients and may thereby help to optimize future clinical decision-making. Our functional analyses suggest the possible involvement of ALCAM in ENZA resistance.

LINC01140 inhibits nonsmall cell lung cancer progression and cisplatin resistance through the miR-4742-5p/TACC1 axis.

Recent studies show that lncRNAs participate in drug resistance and nonsmall cell lung cancer (NSCLC) progression. This study aimed to study the roles and mechanisms of long intergenic nonprotein coding RNA 01140 (LINC01140) in regulating NSCLC progression and drug resistance. Real-time quantitative polymerase chain reaction and western blot analysis were used to detect LINC01140, miR-4742-5p, and transforming acidic coiled-coil 1 (TACC1) expression in NSCLC cells. The interaction between two molecules was examined by luciferase reporter and/or RNA immunoprecipitation assays. Cell invasion, apoptosis, and cisplatin cytotoxicity were assessed by transwell invasion assay, flow cytometry analysis, and CCK-8 assay, respectively. LINC01140 was downregulated and miR-4742-5p was upregulated in NSCLC. LINC01140 inhibited miR-4742-5p expression by competitively binding to miR-4742-5p, while miR-4742-5p targeted TACC1 to inhibit TACC1 expression in NSCLC cells. LINC01140 enrichment repressed the invasive potential and cisplatin resistance and triggered apoptosis, which was reversed by miR-4742-5p overexpression. miR-4742-5p inhibition suppressed cell invasion and cisplatin resistance and accelerated apoptosis in NSCLC cells, while TACC1 silencing abolished these effects. Mechanistically, LINC01140 positively regulated TACC1 expression by sponging miR-4742-5p. In conclusion, LINC01140 inhibited NSCLC progression and cisplatin resistance via functioning as a ceRNA for miR-4742-5p to modulate TACC1.

Regulation of stem cell identity by miR-200a during spinal cord regeneration.

Axolotls are an important model organism for multiple types of regeneration, including functional spinal cord regeneration. Remarkably, axolotls can repair their spinal cord after a small lesion injury and can also regenerate their entire tail following amputation. Several classical signaling pathways that are used during development are reactivated during regeneration, but how this is regulated remains a mystery. We have previously identified miR-200a as a key factor that promotes successful spinal cord regeneration. Here, using RNA-seq analysis, we discovered that the inhibition of miR-200a results in an upregulation of the classical mesodermal marker brachyury in spinal cord cells after injury. However, these cells still express the neural stem cell marker sox2. In vivo cell tracking allowed us to determine that these cells can give rise to cells of both the neural and mesoderm lineage. Additionally, we found that miR-200a can directly regulate brachyury via a seed sequence in the 3'UTR of the gene. Our data indicate that miR-200a represses mesodermal cell fate after a small lesion injury in the spinal cord when only glial cells and neurons need to be replaced.

Significant co-expression of putative cancer stem cell markers, EpCAM and CD166, correlates with tumor stage and invasive behavior in colorectal cancer.

BACKGROUND: The crucial oncogenic role of cancer stem cells (CSCs) in tumor maintenance, progression, drug resistance, and relapse has been clarified in different cancers, particularly in colorectal cancer (CRC). The current study was conducted to evaluate the co-expression pattern and clinical significance of epithelial cell adhesion molecules (EpCAM) and activated leukocyte cell adhesion (CD166 or ALCAM) in CRC patients. METHODS: This study was carried out on 458 paraffin-embedded CRC specimens by immunohistochemistry on tissue microarray (TMA) slides. RESULTS: Elevated expression of EpCAM and CD166 was observed in 61.5% (246/427) and 40.5% (164/405) of CRC cases. Our analysis showed a significant positive association of EpCAM expression with tumor size (P = 0.02), tumor stage (P = 0.007), tumor differentiate (P = 0.005), vascular (P = 0.01), neural (P = 0.01), and lymph node (P = 0.001) invasion. There were no significant differences between CD166 expression and clinicopathological parameters. Moreover, the combined analysis demonstrated a reciprocal significant correlation between EpCAM and CD166 expression (P = 0.02). Interestingly, there was a significant positive correlation between EpCAM/CD166 phenotypes expression and tumor stage (P = 0.03), tumor differentiation (P = 0.05), neural, and lymph node invasion (P =0.01). CONCLUSIONS: The significant correlation of EpCAM and CD166 expression and their association with tumor progression and aggressive behavior is the reason for the suggestion of these two CSC markers as promising targets to promote novel effective targeted-therapy strategies for cancer treatment in the present study.

The novel finding of an FGFR1::TACC1 fusion in an undifferentiated spindle cell sarcoma of soft tissue with aggressive clinical course.

We report an aggressive soft tissue sarcoma with FGFR1-TACC1 fusion occurring in the thigh of an 83-year-old man. Microscopically, the tumor was composed of monomorphic spindle cells arranged in intersecting fascicles. The tumor cells showed ovoid nuclei, fine chromatin, indistinct nucleoli, and elongated eosinophilic cytoplasm. Focal increase in nuclear atypia was noted. Immunohistochemistry showed only focal rare positivity to S100, but negative to SOX10, CD34, STAT6, TLE-1, SMA, and other myogenic markers. An extensive panel of immunostains did not reveal a definite lineage of cellular differentiation. The fusion junction of the chimeric transcript involved FGFR1 exon 16 and TACC1 exon 7, which was similar to those reported in other types of neoplasms such as glioblastomas and epithelial cancers. The transcript was predicted to be in-frame and confirmed by Sanger sequencing after reverse transcriptase-polymerase chain reaction. The patient was treated by marginal excision of tumor without receiving adjuvant therapy. He experienced rapid tumor recurrence with distant metastases and succumbed at 3.5 months after presentation. The finding of FGFR1-TACC1 fusion in a high-grade, undifferentiated spindle cell sarcoma of soft tissue is novel and may offer a potential therapeutic target in the near future.

Xbp1 and Brachyury establish an evolutionarily conserved subcircuit of the notochord gene regulatory network.

Gene regulatory networks coordinate the formation of organs and structures that compose the evolving body plans of different organisms. We are using a simple chordate model, the Ciona embryo, to investigate the essential gene regulatory network that orchestrates morphogenesis of the notochord, a structure necessary for the proper development of all chordate embryos. Although numerous transcription factors expressed in the notochord have been identified in different chordates, several of them remain to be positioned within a regulatory framework. Here, we focus on Xbp1, a transcription factor expressed during notochord formation in Ciona and other chordates. Through the identification of Xbp1-downstream notochord genes in Ciona, we found evidence of the early co-option of genes involved in the unfolded protein response to the notochord developmental program. We report the regulatory interplay between Xbp1 and Brachyury, and by extending these results to Xenopus, we show that Brachyury and Xbp1 form a cross-regulatory subcircuit of the notochord gene regulatory network that has been consolidated during chordate evolution.

Association between TBXT rs2305089 polymorphism and chordoma in Iranian patients identified by a developed T-ARMS-PCR assay.

BACKGROUND: Chordoma is a locally aggressive bone tumor with a high capability of recurrence. Because chordoma often occurs at critical locations next to neurovascular structures, there is an urgent need to introduce validated biomarkers. T-box transcription factor T (TBXT; OMIM: 601397) plays an important role in the pathogenesis and survival of chordoma cells. METHODS: Herein, we aimed to show whether rs2305089 polymorphism is correlated with chordoma in the Iranian population. In order to detect rs2305089, tetra-primer amplification refractory mutation system-polymerase chain reaction (T-ARMS-PCR) was used. In total, 19 chordoma patients and 108 normal healthy individuals were recruited and screened using T-ARMS-PCR. The results were subsequently validated by Sanger sequencing. RESULTS: The genotype distributions and allele frequencies were significantly different among the patient and healthy groups (p-value <0.05). The A allele of rs2305089 showed a significant positive association with chordoma risk (p-value <0.05). DNA sequencing verified the T-ARMS-PCR results as well. This study demonstrated the association between TBXT rs2305089 and chordoma in an Iranian population using a simple, accurate, and cost-effective T-ARMS-PCR assay. CONCLUSIONS: Our results were in line with those of previous studies showing that TBXT rs2305089 is associated with chordoma development. We also developed an efficient T-ARMS-PCR assay to determine the genotype of rs2305089.

A 37†kb region upstream of brachyury comprising a notochord enhancer is essential for notochord and tail development.

The node-streak border region comprising notochord progenitor cells (NPCs) at the posterior node and neuro-mesodermal progenitor cells (NMPs) in the adjacent epiblast is the prime organizing center for axial elongation in mouse embryos. The T-box transcription factor brachyury (T) is essential for both formation of the notochord and maintenance of NMPs, and thus is a key regulator of trunk and tail development. The T promoter controlling T expression in NMPs and nascent mesoderm has been characterized in detail; however, control elements for T expression in the notochord have not been identified yet. We have generated a series of deletion alleles by CRISPR/Cas9 genome editing in mESCs, and analyzed their effects in mutant mouse embryos. We identified a 37 kb region upstream of T that is essential for notochord function and tailbud outgrowth. Within that region, we discovered a T-binding enhancer required for notochord cell specification and differentiation. Our data reveal a complex regulatory landscape controlling cell type-specific expression and function of T in NMP/nascent mesoderm and node/notochord, allowing proper trunk and tail development.

Cell-to-cell heterogeneity in Sox2 and Bra expression guides progenitor motility and destiny.

Although cell-to-cell heterogeneity in gene and protein expression within cell populations has been widely documented, we know little about its biological functions. By studying progenitors of the posterior region of bird embryos, we found that expression levels of transcription factors Sox2 and Bra, respectively involved in neural tube (NT) and mesoderm specification, display a high degree of cell-to-cell heterogeneity. By combining forced expression and downregulation approaches with time-lapse imaging, we demonstrate that Sox2-to-Bra ratio guides progenitor's motility and their ability to stay in or exit the progenitor zone to integrate neural or mesodermal tissues. Indeed, high Bra levels confer high motility that pushes cells to join the paraxial mesoderm, while high levels of Sox2 tend to inhibit cell movement forcing cells to integrate the NT. Mathematical modeling captures the importance of cell motility regulation in this process and further suggests that randomness in Sox2/Bra cell-to-cell distribution favors cell rearrangements and tissue shape conservation.

Plating human iPSC lines on micropatterned substrates reveals role for ITGB1 nsSNV in endoderm formation.

Quantitative analysis of human induced pluripotent stem cell (iPSC) lines from healthy donors is a powerful tool for uncovering the relationship between genetic variants and cellular behavior. We previously identified rare, deleterious non-synonymous single nucleotide variants (nsSNVs) in cell adhesion genes that are associated with outlier iPSC phenotypes in the pluripotent state. Here, we generated micropatterned colonies of iPSCs to test whether nsSNVs influence patterning of radially ordered germ layers. Using a custom-built image analysis pipeline, we quantified the differentiation phenotypes of 13 iPSC lines that harbor nsSNVs in genes related to cell adhesion or germ layer development. All iPSC lines differentiated into the three germ layers; however, there was donor-specific variation in germ layer patterning. We identified one line that presented an outlier phenotype of expanded endodermal differentiation, which was associated with a nsSNV in ITGB1. Our study establishes a platform for investigating the impact of nsSNVs on differentiation.